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I. INTRODUCTION 

The Anderson Impurity problem describes the inter- 
action between a localizaed correlated state and a set 
of uncorrelated orbitals. This ideal system has various 
physical realizations, ranging from the "classical" mag- 
netic impurity in a metal case to the Kondo quantum 
dots nanostructures. In the Kondo regime 1 - a correlated 
many-body state is formed, in which a cloud of extended 
hole states screen the impurity magnetic moment. When 
two of these "impurities" are put close enough the Two 
Impurity Anderson (TIA) problem is set. Magnetic cor- 
relations between the impurities arise, mediated by the 
extended states of the host. For nanostructures, these 
magnetic correlations play a central role in the design of 
miniaturized spin-based devices 2 - made of quantum dots 
(QD). For the "classical" magnetic impurity problem, the 
TIA is a fundamental step towards the understanding 
of heavy fermions systems 3 -^. Man-made nanostructures 
that mimics few Anderson impurity arrays have made the 
full understanding of the correlations present in these sys- 
tems a top priority, because of their possible applications 
in quantum information^^. 

For many years the numerical analysis of Ref.f^] has 
been thought to describe the physics of the two impu- 
rity system in the Kondo limit. As a result the general 
belief 2,3 is that magnetic correlations between impurities 
are generated by the Ruderman-Kittel-Kasuya-Yosida 8 
(RKKY) interaction and that the Kondo effect screens 
such correlations. Very recently it appeare d 4 ' 9 ' 10 : 11 theo- 
retical evidence that some very important aspects of the 
problem was missed in the cited analysis. In fact, it has 
been found in Refs. [ 12 ' 13 ' 14 ] that there is a series of two 
impurity Kondo structures that lead to the formation of 
a ferromagnetic singlet as the ground state of the system. 

The first of that structures corresponds to the forma- 
tion of a Kondo-doublet state, a very simple state medi- 
ated by the interchange of one hole between the impuri- 
ties. That interchange is maximized when the impurity 
spins are aligned, determining the ferromagnetic charac- 



ter of this stage and of the final singlet ground state. The 
relevance of these Kondo-doublet states is two handed, 
they are a constructive proof of the Kondo-doublet inter- 
action and also variational ground states of the spin one 
half (S = 1/2) subspaces. The final stage, which corre- 
sponds to the formation of the singlet ground state, can 
be interpreted as a "composite" single-impurity Kondo- 
singlet, in which what is screened by the formation of the 
Kondo-singlet is the doublet spin, not the impurity ones, 
which remain ferromagnetically aligned. This is called a 
super-singlet in Ref.f 1 ^]. 

Most of the "Kondo" correlation energy associated 
with the formation of the TIA singlet is generated by the 
formation of the Kondo-doublets. The Kondo-doublets 
correlation energy has been shown to be higher than the 
RKKY one for most of the parameter space of the Hamil- 
tonian. Here we analyze the behavior of the Kondo dou- 
blet states beyond the Kondo limit. The effect of high 
energy configurations is analyzed. As those configura- 
tions are directly connected to the vertex of the Kondo- 
doublets, they can be considered order zero corrections 
to the doublets. We found that the Kondo doublet inter- 
action is robust, and its coherence energy is incremented 
by the synergy of the doublets with the new considered 
configurations. The interplay between the different co- 
herence channels that appear when the high energy con- 
figurations are considered can be used to manufacture 
TIA systems of predetermined characteristics. 

This paper is organized as follows: In Section |TT] we 
write down the TIA Hamiltonian and a brief description 
of the Kondo-doublet interaction is given. An analysis 
of the symmetry of the doublets is performed. In Sec- 
tion IIIII we evaluate the odd doublet variational wave 
function (VWF) with the added high energy configura- 
tions. In Sections IIVI and fVl we analyze the response of 
the doublets in the different regimes. In Section fVIII we 
summarize our results and a qualitative interpretation of 
the experimental data of Refs.f 2 ^] is given. In the Ap- 
pendix we discuss some aspects of the variational J n /Ns 
expansion method we used to analyze the system. 
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II. HAMILTONIAN AND KONDO LIMIT 
DOUBLETS 

The Two Impurity Anderson Hamiltonian is the sum 
of the band, hybridization [Hy), and impurity (located 
at Rj = ±R/2 over the x-axis) terms 

H = ^2e ka c\ a c ka + v ^(e* kR] d) a c ka + h.c.) 

k<7 jka 

-E d J2 d) a d 3a + uJ2 4i d iAi d n > W 

where the fermion operator c kG (dj a ) acts on the con- 
duction band fc-state (on the impurity (or QD) at Rj) 
and v = V/y/nb is the hybridization divided by the 
square root of the number of band states. Single state 
energies are referred to the Fermi energy. We renor- 
malize the vacuum (denoted by |F)) to be the conduc- 
tion band filled up to the Fermi energy (Ep), and we 
make an electron-hole transformation for band states be- 
low the Fermi level: &t_ = c ka . We use in the text 

ha 

a "ket" notation for the impurity configurations, the 
first symbols indicate the status of the left impurity 
(the one at x = —R/2) and the second ones the sta- 
tus of the impurity on the right, e.g. |0, f) = d^^F), 
\U)=d{,dUF), \lAV = dl/ Rl dUF). 
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FIG. 1: Some configurations that are present in the Kondo- 
doublet and RKKY interactions; their notation and energy. 
The little black (white) dot corresponds to an electron (hole) 
excitation in the band. The gray circles represent the impu- 
rities. 



In the Kondo limit the impurity level is well below 
the Fermi energy and it can not be doubly occupied due 
to the strong Coulomb repulsion (0 -C Ed U). In 
this regime the two relevant parameters of the TIA are 
the effective Kondo coupling J n = p V 2 /Ed (p a being 
the density of band states at the Fermi level) and the 
interdot distance R. The single impurity Kondo energy 
is given by 8k — D exp (— l/2J n ), where D is the half- 
band width of the metal host. 

In the two impurity case, when the impurities are close 
enough, the higher energy stage of the Kondo nucleation 
is generated by the Kondo-doublet interaction. In the 
Kondo hmrfciS this interaction is described by the follow- 



ing variational wave function 

|A, T ) = \A r ) + vJ2 Z 2 (k)(\A nk ) + \A m )) , (2) 

k 

where the vertex state is given by 

|A T ) = (4 T -4 T )|F) = |0,T>-|T,0) , (3) 
and the states 

\A U k)= bl (e- lk - R / 2 |U)+e +4k ' R / 2 ||,T» , (4a) 

l^m) - &Ix (e-* k - R/2 |T,T>+e +lk - R/2 |T,t» , (4b) 

are the configurations that are obtained applying Hy to 
the vertex state. Minimizing analytically, one obtains 
Z 2 (k) = 1/(S (R) + efe) for the amplitude, and Ed — 
—2Ed — S (R) for the energy of the doublet, where the 
odd Kondo-doublet correlation energy is given by 

S o[e) (R) = D exp(-l/(2±C h (S,R))J n ) , (5) 

where the minus sign corresponds to the even Kondo 
doublet, which is the one generated form the symmetric 
(|0, t) + 1 1> 0)) vertex state. The hole coherence factor is 
given by 

C h (6,R) =J2 Mkh) cos (k h .R) / J2 Z M> ( 6 ) 

h h 

and it is, as a function of R, an oscillating decaying func- 
tion ( \Ch(R)\ < 1, of period ~ Aj? (the Fermi wave- 
length), and Ch(0) = 1, Ch(oo) = 0). Its range, there- 
fore, determines the range of the Kondo-doublet interac- 
tion, i.e. how close the impurities must be in order to 
significatively interact through this mechanism. Direct 
"reading" of the Kondo-doublet state Eq. ([2]) makes clear 
that the Kondo-doublet interaction is produced by the 
interchange of one hole between the impurities using the 
vertex state as a virtual bridge. 

Note that this wave function, Eqs.((2|) to (J6|), is, first, 
a strict constructive proof of the Kondo-doublet inter- 
action, and, second, a variational ground state for the 
S = 1/2 subspaces of the Hamiltonian. For this last 
role, it comprises the lowest energy configurations of the 
subspace and the simplest one, and it makes use of the 
strongest interaction present in the Hamiltonian, there- 
fore we expect it to be a good variational ground state 
of the corresponding subspace. Furthermore, it can be 
taken as the first step in a variational J n /Ns expansion. 
In fact, we follow it two more steps in Ref.fi^] in order to 
analyze the Kondo-doublet RKKY interplay. In the fol- 
lowing sections we analyze the effect of two kind of high 
energy configurations that are directly connected with 
the vertex state, i.e. they can be considered order zero 
corrections to the Kondo-doublets. 

It is wort to discuss briefly the dependence of the 
Kondo-doublets in the particular structure of their ver- 
tex states. If one takes just the |0, f) state as the vertex 
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state, and the corresponding subset of the of \A crLaa k) 
configurations that are connected to it by the hybridiza- 
tion, ones obtains for that VWF an energy —2Ed — §k, 
i.e. no interaction between the impurities. This can be 
confirmed by direct inspection of the VWF, which shows 
that it is the product of a Kondo singlet at the left impu- 
rity "times" a spin up electron in the right impurity. This 
is not strange. If one calculates the corrections to the en- 
ergy of the 1 1, |) and | J., j) states, up to fourth order in 
the hybridization V, one finds no R dependent term for 
them. Instead, if one calculates for the (| f , I) ± | I, t)) 
combinations, one finds that the plus combination gains 
an energy S(i?) (it is the S z = component of the fer- 
romagnetic triplet, an odd state), and that the minus 
combination lost the same amount (it is the even antifer- 
romagnetic singlet). S(i?) is half the RKKY energy. 

Therefore, the particular "structure" of the Kondo- 
doublet vertex states, and that of the RKKY, is nothing 
more than what corresponds given the symmetries of the 
system. Note also that these interactions appear for any 
combination other than the bare ones, (0,1) and (1,0), 
but they are maximal for the proper symmetrized ones 
(1, ±1). This point must be taken into account if one at- 
tempts to make approximations directly over the Hamil- 
tonian, as the Schrieffer-Wolf transformation of Ref. [ 9 ] or 
slave boson methods, otherwise the V 2 Kondo-doublct 
interaction will be missed. 

What is not so trivial, however, is that for the RKKY 
this effect translates into a transfer of energy (±£(i?)) 
between the combinations, whereas that for the Kondo- 
doublet what is transferred is the "connectivity" of the 
vertex state, i.e. its ability to interconnect the Z2 con- 
figurations, the (2 ± Ch{R)) factor in the exponential of 
their correlation energy. This is a consequence of the non- 
perturbative nature of the Kondo-doublet interaction. 

The second stage^^ of the Kondo nucleation corre- 
sponds to the formation of a "composite" Kondo-singlet, 
in which the doublet states play the role that the impu- 
rity states play in the single impurity Kondo-singlet. The 
correlation energy gained in this last stage is given by: 



jk^D exp (-1/(1 TC h (rf,R))J n ) 



(7) 



Thus, the total correlation energy gained by the forma- 
tion of the Kondo singlet in the two impurity case is given 
by E K (R) = &d+1k > 2S K , where S D = max (c^, <5 ) 
is the correlation energy of the dominant doublet. The 
following relation holds Sd{R) > 6k > Jk(R)- For a 
relatively strong "overlap" of the impurities, i.e. 1 > 
\C h (R)\ > 0, one has S D (R) > S K > ^ K (R) and thus 
the Kondo-doublet states are formed at a much higher 
temperature than the super-singlet. 

An extreme case is R = 0, for which Ch = 1- At this 
point the odd combinations of the band states are de- 
coupled from the impurities, and the odd Kondo-doublct 
is the ground state of the system. In the opposite limit, 
R 1— > 00, Ch = 0, one has So = 7if = &k, and the "su- 
per" singlet is just the product of two single impurity 
Kondo singlets, one at each impurity. 



III. KONDO DOUBLETS BEYOND THE 
KONDO LIMIT 

The configurations that generates the Kondo-doublet 
interaction, as analyzed in the previous section, are the 
lower energy ones of the configurations depicted in Figj2l 
i.e. the ones with one electron in each QD (impurity) and 
one hole in the band (plus the vertex state configuration) . 
In the resonant paths the band excitation generated when 
the population of one QD is modified can be reabsorbed 
in any of the two QDs. 
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FIG. 2: Kondo-doublet paths, i.e. configurations connected 
via the hybridization term of the Hamiltonian. The non- 
resonant paths are relevant to determine the full Kondo- 
doublet energy due to the synergy of Kondo structures. 

The band excitation in the non-resonant paths, in- 
stead, can only be reabsorbed at the QD at which it 
was generated. Nevertheless, these configurations must 
be included in the interaction analysis because of the 
strong synergy of the non-perturbative Kondo structures 
(Refs.[ 8 ], pag.155 ,f^],pag.53). In the following we ana- 
lyze the effects of the high energy configurations shown 
in FigfJl the ones with no electrons in the QD's and the 
ones with a double occupied QD, on the Kondo-doublet 
states. These configurations are directly connected to the 
vertex state. 

Therefore, rewriting Eqs.(|3"l fla)) for completeness, the 
properly symmetrized one-electron in-two-QDs configu- 
ration that acts as vertex of the odd Kondo-doublet is 



|Ar> = (4a-Ol^) = M-kT,0) , 



(8) 



of energy —Ed- This configuration acts as a bridge be- 
tween the other components of the doublet, and it is a 
"virtual" state in the Kondo limit. Particularizing to the 
spin-up component of the doublet, the hybridization Hy 
connect the |^4|) configuration with the following ones 
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of energy —2Ed+ek, by promoting an electron from below 
hp to the empty QD. In Refs. [ 12 i 13 ' 14 ] just these configu- 
rations were considered. The first one is a non-resonant 
configuration and the second one is the one that gen- 
erates the strong ferromagnetic correlation between the 
impurities. This can be checked by the evaluation of their 
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connectivity factors, i.e. the configuration to vertex state 
matrix element of the hybridization, normalized by — 2v 

(A^HvlAnk) ~ 1 , (10a) 
(A-\\Hv\A^k) ~ (1 + cosk.R), (10b) 

this last element depends on the inter-QD distance i?, 
it generates correlations between the status of the two 
QDs. 

The hybridization also connect the \A-\) state with the 
following ones: the resonant 

k T > = (-4 e^ q - R/2 \F) + cJ T e + ^- R / 2 |F)) (11) 
configuration, of energy e g and connectivity 

(A t |iJ v |g T ) ~ (1-cosq.R) , (12) 
and the non-resonant 

\U kr ) = &I T (e- ik - R/2 |0, IT) - e +4k ' R/2 | IT, 0)) (13) 
configuration, of energy —2E d + U + e k and connectivity 

(Ai\H v \U m ) ~ 1 . (14) 

Including all these states, the variational wave function 
for the study of the odd doublet is given by 

\D o] ) = L4 T )+v^ Z 2 (k)(\A nk ) + \A m )) 

k 

+ v]T Z (q)\ qi )+vJ2 Zu(k)\U k1 ) , (15) 

q k 

where the k (q) sums are over hole (electron) excitations. 
In the following we will use the variational amplitude fac- 
tors to refer to the corresponding set of configurations. 
Minimizing the expectation value of the odd doublet en- 
ergy, we obtain 

Z 2 {k) = -l/{E + 2E d -e k ) , (16a) 
Mo) = -l/(Eo-e q ) , (16b) 
Z v {k) = -l/(E + 2E d -U-e k ) , (16c) 

and 

E = -E d + v 2 J^(2 + cos k x R)Z 2 (k) 

k 

W 5^(1- cos q x R)Z (q) + v 2 £ Z v {k) , (17) 

q k 

where E a is the energy of the odd Kondo-doublet. As- 
suming E Q = —2E d — 5 0l we rewrite the amplitude factors 
as 

Z 2 (k) = l/(S + e k ) , (18a) 
Z (q) = \/{2E d + 6 + e q ) , (18b) 
Zu(k) = l/(U + 5 + e k ) , (18c) 



and the equation for the doublet energy becomes a self- 
consistent equation for the doublet correlation energy 6 



E d + S 
V 2 Po 



I h {S )(2 + C h (S )) + I h (U + S ) 
+ I e (2E d + 5 )(l - C e (2E d + 6 )) , (19) 



where the Kondo integrals are given by I x (uj) = 
(1/nbPo) Y] x l/(^ + ex) an< i the hole (electron) coherence 
factors by 



C x {co,R) 



1 . cosk^.R , , , . . . 

(zZ —TZ ) / J xM- (20) 



n h p 



lu + e x 



Further elaboration of Eq. (fl9| depends on the particular 
system and regime in study. The Kondo limit in 1, 2, and 
3 dimensions was analyzed in Ref.fi^], obtaining Eqs.([5l 
[7]). In the following we center our analysis in the finite 
E d , U effects generated by the high energy configurations, 
while using a simple model for the extended k states. 
With the half-filled flat band model, of half-band width 
D, the Kondo integrals result to be 



I e (w) = I h (w) = (a;) 



: hi 



D + oj 



ui 



(21) 



For the coherence factors, further choices must be made. 
A one dimensional system^ has the advantage that deco- 
herence comes only from the energy spread (~ 8 ) of the 
involved excitations and not from the angular terms that 
are characteristic of higher dimensions 12 . In near 2D QDs 
systems made in semiconductor heterostructures such an- 
gular dispersion effects can be minimized by electron- 
focusing technics, rendering the system response close to 
the ID results. Furthermore, in ID the coherence factors 
can be analytically evaluated, which is very convenient 
for the purposes of the present study. Thus, particular- 
izing for a ID band system, the coherence factors are 

C h{e) {u>, R) = [cos (r* ± uj*)[Ci(r* + lu*) - Ci(u*)} 
±sin(r* ±w*)[Si(r* +lo*) - Si(w*)]]/7jr(w) , (22) 

where r* = k F R — 2irR/\ F , to* = k F R lu/D, and Ci 
(Si) is the Coslntegral (Sinlntegral) function. Both Ch 
and C e tend to cos (k F R) for lo i— ► 0, i.e. when the 
band excitations that generate the correlation are the 
ones with k ~ hp. 

In FigJ3]we plot the coherence factors. It can be seen 
that the electron coherence factor has a period a lit- 
tle lower than X F because it involves excitations with 
k > k F , the opposite is true for the hole coherence factor. 
This "mismatch" between the coherence factors opens 
the possibility to design arrays with different character- 
istics, as we will discuss in the next sections. 

For the even Kondo-doublet, which vertex is the sym- 
metric l^) = [d ] Ra + d[ a )\F) = \0,a) + |er,0) state, 
it holds the same equations but with the Q x i— * — Q x 
change. 

The odd doublet variational wave function Eq. lTSl) is 
valid for the analysis of the < 2E d < U region, in 



5 



and the norm 



8 

e 



1.00 
0.75 
0.50 
0.25 
0.00 
-0.25 
-0.50 
-0.75 
-1.00 



A 

M 

A / 


."•^cosCkpR) ; ". 

r — t\ — n — ^ — L 


A >■ 

A ''■ 


/ A ' ■/ V A /; 



0.0 0.5 



1.0 



1.5 



2.0 



2.5 



interdot distance R / X„ 



FIG. 3: Hole and electron coherence factors as a function of 
the inter-impurity distance, for u = 0.001 D. In ID their 
range is given by the Kondo length. 



which an empty impurity state (of energy 0) is the main 
channel for the hybridization to act. For the U < 2Ed < 
2U region, in which the main hybridization channel is 
a double occupied impurity (of energy — 2Ed + U < 0), 
a similar doublet generated by electron-hole symmetry 
arguments can be used. The vertex state of such doublet 
is given by 



I^t) = IIT,T)-|T,IT) , 



(23) 



and the corresponding equations and properties are sim- 
ilar to the ones analyzed here. 



(D\D) = 2 [ 1 + v 2 { E (2 + cos k.R) Z 2 (k) 2 



E z ^ 2 



E (l-cosq.R) Z (q) 2 }] , 



k q 

of the doublet state. They appear two new sums, 

1 

n bPo ^ (w + e x ) 2 



J x (u>) 



and 



D x (u,R) = ( 



nbPo 



E 



cos k x .R 

(u> + e x )' 



(26) 



(27) 



(28) 



that can be easily evaluated as derivatives with respect 
to ijj of the previously defined I x and C x . With the as- 
sumed band model it results Jk{u) = D/lj(ui + D). For 
D x (lj, R) it results a large expression, suffices to said that 
its behavior is very similar to that of C x (ui,R); its am- 
plitude decays a little slowly because its weighting factor 
(Z 2 ) is more biased towards the k ~ ftp excitations than 
the one of C x (Z x ). Thus the spin-spin correlation is 
given by 



(S L .S R ) = (D\S L .S R \D)/(D\D) = 

3 V 2 Po D 2 J 2 

4 + 4 V 2 Po (J 2 (2 + D 2 ) + J V + J (l - D )) 



(29) 



where J 2 = J h (S ), J v = J h (U + S a ), and J = J e (2E d + 
S ) and the same convention is used for the D x (uj) argu- 
ments. 



A. Impurity Spin-Spin correlation 

To characterize the response of the TIA system in the 
different regimes we evaluate some observables. The most 
relevant is the spin-spin impurity correlation (Sl-Sr), 
which clearly signals the extended nature of the Kondo 
doublets. The mean population of one impurity indicates 
the changes from Kondo to intermediate valence (or An- 
derson) regime. At this effect we need to evaluate the 
following quantities: the expectation value of the spin- 
spin correlation, 

(D\S L .S R \D) = I v 2 E cosk - R Z 2(k) 2 , (24) 
k 

of the average population of one impurity 

(D\n R \D)= (D\n Br +nni\D)= 1 + 2 v 2 
{E(2 + cosfc.i?)Z 2 (£0 2 + E Z v (k) 2 }, (25) 

k k 



IV. U EFFECTS. 

Corrections due to a finite value of U are the lower 
ones in the Kondo doublet VWF (Eq. (fr5"|) 1, for its range 
of applicability. For < 2Ed < U, it results that 
Z 2 ^> Zq ^> Zjj . As these corrections come from the 
non-resonant Zjj configurations they have little depen- 
dence on R. Working over Eq. fTT?]) . and considering only 
the Z 2 and the Zjj configurations effects, one obtains 

6 (U)/6 = (1 + ^/[/) 1 /(2+c,(5 )) i (30) 

which gives the incremental ratio of the Kondo doublet 
energy. It depends very weakly on J n and R, through its 
dependence in Ch- The incremental band is plotted in 
Fig|4]as a function of D/U. 

A remarkable aspect of this correction is its use of the 
synergy of Kondo structures, i.e. the strong increment 
of the doublet coherence energy related to the increment 
of the number of available states for the jumping exci- 
tation. This effect is achieved with little effects in other 
properties of the doublet. At D/U ~ 0.2 roughly a ten 
percent increment of the coherence energy is obtained, 
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FIG. 4: 5 (U)/5o as a function of D/U. The relative change 
is greater the lower is Ch, covering the gray zone in the figure, 
for C h = 1 to 0. 



whereas that the relative weight of the Zu configurations 
in the VWF, given by E^/2£Zf = Jk(U)/2J h (S ) 
~ 5 D/2U 2 , is of the order of 1CT 4 . Thus the spin- 
spin correlation generated by the Kondo doublet remains 
nearly unchanged. 



V. E d EFFECTS. 



i.e. the same equation than for the single impurity prob- 
lem, Ref.fi, Eq.(7.17)] but with a connectivity factor of 
3. For the even doublet it is obtained 



E d + 6 e = V 2 p Q {ln- 



D 2E d + S e 
h 2 in ■ 



D 



2E d 



) • (32) 



See that the Zq configurations are decoupled from the 
odd doublet (at R = 0, C e = 1), whereas that they make 
an important contribution to the even doublet, which in 
turn has its connectivity to the Z2 configurations reduced 
to 1. 
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FIG. 5: Odd doublet coherence energy, Eq. (|19[) . as a function 
of E d for 7? = 0.0 and J n ,J D = 0.0732. For E d > D (Kondo 
regime) we keep constant the V 2 /Ed ratio, whereas that for 
D > Ed (Anderson regime) the V 2 /D ratio is maintained 
constant. 



We separate the study of Ed corrections on the Kondo 
doublet properties in two regimes. For Ed > D (Kondo 
regime) we use the usual effective Kondo-coupling ap- 
proximation, keeping constant the V 2 /E d ratio, i.e. we 
study our equations for fixed values of J n = V 2 p /Ed- 
For D > Ed > (intermediate valence or Anderson 
regime) the hybridization amplitude (V) is keep con- 
stant, thus in this region we study the system response 
as a function of Jd = V 2 p /D. Due to the structure of 
Eq. lfTT))) . the one that determines S for given Ed, U, V, R 
values, the U effect analyzed in the previous section is 
nearly the same multiplicative factor whatever the value 
of Ed, thus we disregard the Zu configurations in this 
section. 

One of the effects that arise in considering a finite value 
for Ed is that involved in the usual Ed, fl » J approx- 
imation that allows for the "exponential" solutions of 
Eq. (fl"9]) . Eqs.([5]). To visualize this point we first solve 
Eq. (fT!))) as a function of Ed at R = 0. Although this can 
be seen as an unrealistic situation for the "classical" mag- 
netic impurity problem, it is a very important setup for 
quantum dots laterally coupled to a quantum wir o 6 ' 11 . 
At this value of R one has Ch,C e = 1, thus the equation 
for the coherence energy of the dominant odd doublet at 
R = results 

E d + 6 o = 3V 2 p ]n^±^ , (31) 



In FigJ5]we plot the odd and even doublets coherence 
energy as a function of Ed , for R = and J n , Jd = 
0.072382. This value of J n corresponds to a Kondo en- 
ergy equal to 5k = 0.001 D and, in the Kondo limit, to 
an odd doublet energy of 5 3 =S o (R = 0) = 0.01 D. The 
corresponding even doublet energy is 5i — S e (R — 0) = 
0.000 001 D. Lowering Ed from the Kondo limit it can 
be seen that 8 decreases a little from his "exponential" 
expression S = D exp (— 1/3 J n ), but this is still a good 
approximation up to very low values of Ed/D, see the 
point line in the left panel of Fig0 This point line gives 
the solutions of Eq. (f3"Tj) maintaining V 2 /Ed constant in 
the Ed/D < 1 region. The full line in the left panel 
is for a constant value of V, such that Jd — J n Ed/D 
remains constant, which is more in line with the exper- 
imental situation for quantum dots. The dashed line is 
the even doublet coherence energy which is equal to that 
of the odd doublet for Ed — 0, i.e. when the Zq configu- 
rations have the same energy than the Zi configurations. 
Increasing Ed makes the Zq configurations a more expen- 
sive path, and thus 5 e quickly decreases. 

In FigJS]we plot the mean population of one impurity 
for the same situations depicted in Figj5] These values 
are very similar to the ones corresponding to the sin- 
gle impurity casei. Note that the "Kondo coupling ap- 
proximation" , point line, gives higher occupation values 
than the "constant-^" curve. This is because the former 
corresponds to lower values of the hybridization, such 
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FIG. 6: Mean impurity occupation as a function of Ed, for 
the same parameter values than the previous figure. The main 
characteristic of the Kondo ((rid) — 1) and Anderson (or in- 
termediate valence) ((rid) < 1) regimes can clearly be seen. 



that V 2 /Ed remains constant. For the even doublet, as 
measure that Ed decreases, the mean occupation remains 
high (~ 1) because of the reduced connectivity of its Z 2 
configurations until the weight of the Zq configurations 
becomes relevant, and then the mean impurity occupa- 
tion of the even doublet becomes lower than the one of 
the odd doublet. 



A. Kondo regime 

In this region (Ed > D) the usual Ed,D^> 5 approxi- 
mation is a good one, and we use it to recast the solutions 
of Eq. (|3"T|) in an "exponential" form. Depending on the 
value of J n (i.e. V) this is a good approximation down 
to even lower values of Ed/D, see the point lines in the 
left panel of Figs|5][Hl In this regime the effect of the Zq 
configurations in the odd doublet is similar to that of the 
previously discussed Zjj configurations. For the relative 
increment of the Kondo doublet energy one obtains 

S (E d )/S = (1 + D /2Ed) {1 - c ^ 2Ed))/{2+Ch{5 ^ , (33) 

the main difference with the Zjj configurations case is 
that the connectivity of the Zq configurations in the odd 
doublet depends on the electron coherence factor. This 
is reflected in the numerator of the exponent in Eq.(33), 
I — C e (2Ed), which reduce the size of these corrections 
compared with those induced by the Z\j configurations. 
Again, this increment in the coherence energy of the dou- 
blet is obtained with a very little participation of the in- 
volved configurations in the total weight of the doublet 
wave function, and thus other properties as the mean im- 
purity population or the impurity Spin-Spin correlation 
remain nearly untouched. 
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FIG. 7: (Sl-Sr) as a function of Ed, for the same parame- 
ters than the previous figures. In the Anderson regime the 
ferromagnetic response of the odd doublet is reduced by the 
increasing weight of the vertex state. For the even doublet it 
is also important the contribution of the Zq configurations. 



In Fig|7]we plot the Spin-Spin impurity correlation for 
the same situations depicted in Figs 051 This, of course, 
is a characteristic of the two impurity Kondo doublets 
that has not analogy in the single impurity case. The 
ferromagnetic response of the dominant odd doublet re- 
mains close to its saturation value (~ 1 /4) up to low val- 
ues of Ed/D, where the weight of the vertex state starts 
to be significative. For the even doublet the same effects 
that determine its mean impurity population curve can 
be traced in its Spin-Spin correlation. 



B. Anderson regime 

In this region (Ed < D) there is a rich new structure 
in the Kondo doublets behavior due to the influence of 
the Zq configurations and their electron driven coherence 
factor. The relative weight of these configurations in the 
doublet processes increases as measure that Ed decreases, 
i.e. when their energy becomes comparable to that of 
the Z2 configurations, making them a not too expensive 
path for the jumping excitation. Moreover, these paths 
are resonant, their connectivity depending on the elec- 
tron coherence factor. As the "periods" of the electron 
and hole coherence factor are respectively a little shorter 
and a little longer than Xp, a fine tuning of the parame- 
ters of the system, in the quantum dots heterostructures, 
allows for the manufacture of arrays with different char- 
acteristics. 

In Fig(5] we show the Kondo doublets coherence en- 
ergy, mean impurity population, and Spin-Spin correla- 
tion for a value of Ed in the intermediate valence region, 
E d = 0.2 D, and J D = 0.0723824. Due to the influ- 
ence of the Zq configurations the maximum of the odd 
doublet energy is not at R = 0, but at R ~ 0.25 Xp, a 
distance that in the Kondo limit (Ed 3> D) is close to the 
first odd-even doublet transition. In the Kondo limit the 
coherence energy of the odd doublet decrease, as mea- 
sure that R increase, because the connectivity of the Z2 
states, given by (2 + Ch(S, R)), decreases. In the present 
situation this fact is overcompensated by the increase of 
the Z states connectivity (I — C e (2Ed + S,R)), given 
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FIG. 8: Coherence energy, mean impurity population, and 
Spin-Spin correlation of the Kondo doublets as a function 
of the inter-impurity distance R, for Ed ~ 0.2 D, and 
Jd = 0.0723824. The dashed vertical lines mark the odd- 
even doublet transitions. 



FIG. 9: Coherence energy, mean impurity population, and 
Spin-Spin correlation of the Kondo doublets as a function of 
the inter-impurity distance R, for Ed — 0.5 D and Jd = 
0.05428. At this value of Ed the response of the system is 
already similar to that of the Kondo limit. 



that C e oscillates quickly than Ch- Thus the Zq configu- 
rations, decoupled from de odd doublet at R = 0, sustain 
a longer first dominated region for the odd doublet and a 
strong maximum at R ~ 0.25 A^. These configurations 
also generates strong oscillations in the mean population 
of the impurities, because they have a null impurity pop- 
ulation, see Fig[2] And, given that the Spin-Spin correla- 
tion depends mainly in Dh(— Ch), Eq. (|29|) . the extended 
range of the odd doublet makes it to include a slightly 
antifcrromagnetic region near the first odd-even doublet 
transition. 

In Figl5] we show the response of the system for a 
greater value of Ed than in previous case. For Ed 
0.5 D, and Jd = 0.05428, the general behavior of the 
system, as a function of R, is already similar to that of 
the Kondo limitAS. The mean population of the impuri- 
ties is closer to one and has lower amplitude oscillations 
than in the previously analyzed case. The value of the 
ferromagnetic correlation of the dominant doublet is also 
closer to its saturation limit. The effect of the Zq states 



can be seen in the still extended first odd doublet region. 

VI. CORRECTIONS IN THE RKKY 

The inclusion of high energy configurations in the anal- 
ysis of the interactions present in the Two Impurity An- 
derson Hamiltonian also generates fourth order terms 
others than the usually discussed RKKY interaction. In 
FigHOlwe depicted all the paths that generate fourth or- 
der correlations between two aligned Anderson impuri- 
ties. The dashed curved lines mark the RKKY path. The 
upper most path, that goes through the full empty impu- 
rity configurations, is the superexchang o 15 i 16 (or double 
exchange) one. The lower ones, below the horizontal line 
at the middle of the figure, are the "hole-electron" sym- 
metric of the upper ones. 

Although usually ignored in applications of the RKKY, 
these terms has been known for a long time. They were 
reported in Ref.£i£], and carefully analyzed in Ref.fi&J. 
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FIG. 10: Paths corresponding to the FM RKKY-like inter- 
actions. Fourth order paths that only contribute to single- 
impurity corrections are not shown. The dashed lines mark 
the RKKY path. 



The main result is that the superexchange usually com- 
petes with the RKKY, i.e. they have opposite signs, and 
eventually it dominates the RKKY for short distances be- 
tween the impurities. These perturbative interactions are 
additive, there is not synergetic effects for them. More 
recently, they have been re obtained in Ref. fa]. 



VII. CONCLUSIONS 

We analyze the effect of high energy configurations on 
the Kondo doublet interaction. We find analytic expres- 
sions for the corrections to the Kondo doublet coherence 
energy around the Kondo limit. We show that in that 
limit a relatively strong increment of the energy is pro- 
duced without significatively affect others properties of 
the doublets. These corrections, due to the full empty 
and double occupied impurity configurations, are multi- 
plicative. 

In the intermediate valence region we find that, as 
measure that full empty impurity configurations becomes 
energetically accessible, there is an interesting interplay 
between hole and electron driven coherence effects. The 
maximum of the odd doublet is shifted from R — to- 
wards R ~ Xp/4 and the first R- region dominated by the 
odd doublet is significatively increased, making a slightly 
antiferromagnetic region accessible. 

As pointed out in Ref. fa], the experimental results of 
Craig et aia can be interpreted as the resonances cor- 
responding to the odd and even Kondo doublet states. 
But accordingly with our previous results^, we can not 
rule out that the lower resonance corresponds instead to 
the formation of the two impurity Kondo super-singlet. 
Our results for the odd-even Kondo doublet transitions 
as a function of the inter-impurity distance are compat- 
ible with the data of Wahl et aU> for Co adatoms on a 



cupper surface, a detailed analysis will be present else- 
where. We must pointed out that both the authors of 
Ref. fa] and Ref. fa] interpreted their results in terms of 
the RKKY exchange interaction. 

Therefore, we have shown that the Kondo doublet in- 
teraction is robust beyond the Kondo limit. The inter- 
play between the different coherence channels can be used 
to manufacture two impurity systems of predetermined 
characteristics. 
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APPENDIX: ON THE J n /N s EXPANSION 

To analyze the Two Impurity Anderson system we use 
an extension of the "variational 1/Ns expansion", Ref. f 1 ] 
pag.223 (hereafter HW), to the few impurity cas e 12 ' 16 . 
This method is based in the Kondo singlet VWF designed 
by Varma and Yafet^a. It was later used by Gunnarsson 
and Schonhammcr— to successfully explain the spectra 
of rare earth compounds, and has become known as the 
Variational 1/Ns Expansion, where N$ is the degeneracy 
of the localized orbital. Although detailed in Hewson's 
book, we point here to some aspects of the method not 
covered there. 

First, the connection of the method with the Varma 
and Yafet singlet, see HW-Eq.(8.44), that is the VF sin- 
glet. VF work is not cited in HW. 

Second, the expansion parameter is J n /Ns, not just 
1/Ns- See HW-Fig. 8.20. Each new generation of states 
that are included in the variational wavefunction is com- 
prised of the new configurations that are obtained from 
the previous one by the application of Zfy. Therefore, 
calling generation zero to the starting \F) state (taken 
also as the energy reference), odd generations have one 
electron in the impurity and their energies start from 
—Ed , and even generations have zero (or two) electrons 
in the impurity, and their energies start from zero (or 
—2Ed + U > 0). In the Kondo limit the even generations 
play the role of nearly virtual states connecting the odd 
generations states (their total weight in the VWF is much 
lower than that of the odd generations). Thus, a factor 
p V 2 /Ed = J n appears every two generations. Typical 
values of J n are 0.1, for 8k — 0.01D, to 0.05, for which 
8k — 0.0001-D. The results of this method are exact in 
the J n /N s -> limit. 

Third, range of application. Taking into account the 
VF VWF, and the first correction to it, the approximate 
ground state for the single impurity Anderson system is 

ka kqa 
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the variational amplitudes result to be Z k = v/(5K + ek), 
and Y kq = v Z k /(—E s + e k + e q ). The energy of the 
Kondo singlet is given by 

E s = -E d - 5 K + Ej , (A.2) 

where the Kondo energy 8k comes from the resonance 
between the Z k configurations, which use the Fermi sea 
configuration as a nearly virtual bridge. Ej is the single- 
body impurity correction, equal to — J/2 for the level of 
approximation used in Eq. ljA.ip . It comes from the in- 
terplay between each Z ko configuration and their derived 
Y ko q configurations. Therefore Eq. (|A.2j) ilustrates the 
versatility of this variational method: it gives both the 
non-perturbative 8 Kondo term as well as the standard 
perturbative ones. At this level, if one wants to compute 
corrections in U , one must also include the configurations 
with two electrons in the impurity, &Lbp-f| T4)i t na t, are 
obtained from the Z k states, these configurations are not 
shown in HW-Fig. 8.20. In this way the VWF can be 
expanded to the desired level of accuracy, depending on 
the situation to analyze, see for example the works by 
Gunnarsson and Schonhammer. If few generations are 
needed, results can be obtained in a nearly analytical 



way. In other cases, the variational amplitudes can be 
determined by numerical methods. 

The picture above is for the "usual" Kondo regime 
(U ^> Ed ^S> 0), in which the main "virtual" configura- 
tions are the ones with zero electron in the impurity, and 
the ones with two electrons are relegated to a secondary 
role because of their energy cost (—2E d + U 0). For 
the Ed < U < 2Ed region, in which the single occu- 
pied impurity is still the lower energy configuration, but 
the double occupied one is the main "intermediate" state 
(—Ed < —2Ed + U < 0), an electron <-> hole mapping 
must be done. In this region the VF singlet reads 

\S u K ) = \][)+Y,Z q c\ a \a) , (A.3) 

qa 

and it can be expanded to any desired level of accuracy 
in the same way than in the previously discussed region. 

Therefore, this method can be applied almost anywhere 
in the space of parameters. The U = 2Ed can not be 
directly accessed, but it can be approached from above 
and below. This method can be easily extended to the 
few impurities case o 12 i 16 . 
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